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We study the cavity field's and atomic asymptotic mean excitation numbers due to anti-rotating 
term (AMENDART) in the circuit Quantum Electrodynamics (circuit QED) system, composed 
of a two-level atom and a single cavity field mode, subject to Markovian damping and dephasing 
mechanisms. We show that the AMENDART are above the thermal values, and their behavior 
is analyzed analytically and numerically for typical parameters in circuit QED implementations 
described by the Rabi Hamiltonian. We point out that "parasitic elements" such as other cavity 
modes or eventual off-resonant atoms, also contribute substantially to AMENDART. 
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I. INTRODUCTION 

Cavity Quantum Electrodynamics (cavity QED) is one 
of the most active areas of research in quantum optics 
and quantum information nowadays. It deals with the 
light-matter interaction for electromagnetic (EM) fields 
with a few photons and a small number of real or artifi- 
cial atoms confined in optical and microwave resonators. 
It can be implemented in a variety of architectures, such 
as atomic ensembles, Rydberg atoms, quantum dots, su- 
perconducting circuits, Bose-Einstein condensates, polar 
molecules, trapped ions, etc [3-|B|. Here we focus on the 
cavity QED realizations in superconducting circuits, the 
area known as circuit QED 0, H[ , where the one- photon- 
one-atom interaction is currently realized with different 
types of artificial atoms, and interaction between differ- 
ent atoms, placed deterministically inside the cavity, is 
realized by means of the cavity field acting as the quan- 
tum bus. Some advantages of using superconducting cir- 
cuits are the possibility of engineering the properties of 
the artificial atoms, such as the energy levels configu- 
ration and the atom-field coupling, rapid tuning of the 
cavity's and atom's frequencies and the freedom in fabri- 
cating several artificial atoms at specific locations within 
the cavity. Besides studying fundamental problems in 
quantum mechanics, such as entanglement, dissipation, 
decoherence and the measurement problem, circuit QED 
is also an important tool for testing basic quantum al- 
gorithms M, eng ineering nonclassical states of light and 
matter [llul, llOj ]. achieving the strong and ultra-strong 
light-matter coupling 1 1 lL 1 1 211 . implementing rapid time- 
dependent phenomena [la - l2C| and novel technologies for 
quantum information processing [H, Hi]], etc. 

The simplest Hamiltonian, deduced from first princi- 
ples, that describes the interaction between a two-level 
atom and a single mode of the quantized electromagnetic 
field in a cavity is the Rabi Hamiltonian (RH) [22j]. It 
reads (h = I) 

H R = n + QE + gpo-y, (1) 

where we set the cavity frequency to v = 1, the atomic 
transition frequency is f2 and g is the atom-field cou- 
pling constant. The cavity field quadrature operators 



are x = (a^ + a) /y/2 and p = i(a j( — a)/ \f2 1 a and are 
the bosonic annihilation and creation operators, satisfy- 
ing the commutation relation [a,a^] = 1, and n = a^a 
is the photon number operator. The atomic operators 
are E = \e) (e|, a x = cr_ + <r + , o~ y = i(a_ — a + ) 
and o- z = E- | g) (g |, where er + = |e) (g\ and ct_ = n\ , \g) 
and |e) denoting the ground and excited atomic states, re- 
spectively. The exact diagonalization of RH has not been 
achieved yet, so to obtain analytical results one usually 
performs the Rotating Wave Approximation (RWA) by 
neglecting the anti-rotating term (ART) (aer_ + a^o- + ), 
responsible for simultaneous creation/annihilation of one 
photon and one atomic excitation [2^, [24| ■ In this case 
one gets the celebrated Jaynes-Cummings Hamiltonian 
(JCH) [2 5L Hfjj], which has a simple analytical solution 
[1, [H, |24| that led to prediction of several pure quantum 
effects, many of them already verified experimentally [l|- 

i sans Bui. 

The RH or JCH alone do not describe the current sit- 
uation in circuit QED because of dissipation and deco- 
herence arising from the system interaction with cavity's 
and atomic environments. So instead of the Schrddinger 
equation one has to use the quantum master equa- 
tion for the density operator p, whose general form is 




j t P + i[H,p] =Cp, (2) 

where H is some effective Hamiltonian and the superop- 
erator C is the Liouvillian that accounts for the influence 
of the environment. In most situations the Liouvillian 
consists of three parts, C = C K + C\ + £ 7 , where the su- 
peroperator C K (C\) describes the cavity (atom) damping 
by the thermal reservoir with the mean photon number 
n, and n (A) is the cavity (atom) relaxation rate that can 
be determined experimentally 0, H, H, G1, H3, [H, [H| . An- 
other common source of decoherence in superconducting 
circuits is the pure atomic dephasing represented by £ 7 , 
that arises due to low- frequency 1 // noise [Hj], and 7 
denotes the pure dephasing rate (3(| H|[ . The preferred 
theoretical pick to describe the majority of experiments, 
owing to its simplicity and wide applicability range, is 
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the "standard master equation" (SME) H3, IM| 

C K = k{u+ l)V[a] + KfiV[<J] (3) 
C x = X(n + l)V[a-}+XnV[a+] (4) 

A = l^Wz], (5) 

that can be deduced microscopically by making the Born- 
Markovian approximation on the system-reservoir inter- 
action [30], where the so called Lindblad superoperator 

V[$\p = i (2$p$ f - $ f $p - p$ f $) (6) 

preserves the hermiticity, normalization and positivity of 

Within the quantum trajectories approach a Marko- 
vian reservoir represented by a Lindblad kernel can be 
viewed as a "detector" performing continuous measure- 
ments on the system [23 . [34], [35| . Considering, for in- 
stance, £ =x£'[ ( I > ] and rewriting ([2]) as dp/cfa — Lp = Rp, 
where Rp = x^P^ an d Lp = — i[H, p] - {<frt<&, pIy/2, 
one may check that the formal solution for p(i) is [30| 

p(t) = I dt m \ dt m _i ■■■ dtip m (7) 
m=Q Jo Jo Jo 

p m = eHt-t^feHtn-tn-j ... Re Lt lp ( Q y (g) 

This procedure decomposes the quantum dynamics con- 
tained in the master equation into an infinity of quantum 
trajectories, where each occurrence of R corresponds to 
an instantaneous quantum jump of the system's state 
35], and the exponentials exp[L (^fc ~ ifc-i)] describe 
the non-unitary system evolution between the quantum 
jumps. Thus p m can be interpreted as the evolved (non- 
normalized) density operator conditioned to quantum 
jumps at times t%, t<i . . . t m , and its trace gives the proba- 
bility of this particular sequence. Such an interpretation 
describes the Markovian dissipation as continuous mea- 
surement of the system by the reservoir [34j, where the 
quantum jumps describe the "clicks" of the "detector" 
[35j , and the exponential superoperators describe the sys- 
tem evolution without "clicks" but still under continuous 
monitoring |29i. 130}. 

The JCH shows good agreement with experiment and 
with numerical calculations in the majority of problems 
of practical interest [28j |. However, in some occasions the 
ART becomes important, so many studies investigated 
the range of validity of RWA and new physics beyond 
RWA (see references in [36j|). In particular, in [37j was 
shown that the combined action of the anti-rotating term 
and Markovian dissipation leads to incoherent photon 
creation from vacuum, so the asymptotic mean photon 
number is slightly above the thermal value. A detailed 
analysis of this phenomenon was performed in [36j . in 
[2H the cavity emission rate due to ART was calculated, 
and the relation of ART to the entropy operator was de- 
scribed in 39]. Moreover, in the study [40] the errors 
in the zero-excitation state preparation due to ARTs in 



two-atom Markovian cavity QED system were estimated, 
while the creation of transient entanglement between two 
atoms due to ARTs was reported in |4lj . 

Here we study analytically and numerically the cavity 
field's and atomic asymptotic mean excitation numbers 
due to anti-rotating term (AMEND ART) for the current 
parameters in circuit QED, taking into account common 
"parasitic" elements, such as other cavity modes and off- 
resonant two-level atoms, showing that in many cases 
these elements contribute significantly to the generation 
of atomic and cavity field's excitations due to ART. 

II. ANALYTICAL AND NUMERICAL RESULTS 

In this section we analyze analytically and numerically 
the asymptotic mean numbers of photons and atomic ex- 
citations, (n) and (E), respectively, in the circuit QED 
setup described by the master equation @ with H = Hr 
and dissipative kernels ([3])-([5]). Since RH cannot be inte- 
grated exactly, to obtain approximate analytical results 
we perform the fc-photons approximation (denoted below 
by the subindex k) by assuming that at most k photons 
are present in the cavity. The asymptotic values are eval- 
uated by equaling to zero the left-hand-side of the result- 
ing Heisenberg equations of motion. Under the 1-photon 
approximation, valid for (ni) <C 1, for zero-temperature 
damping reservoirs (n = 0) we get simple formulae for 
the AMENDART [36| 

(m) = ^(2G + Xs) (9) 
(Ex) = £(2G + ks), (10) 

where G = g 2 T, T = 7 + (k + A)/2, s = A 2 + T 2 , T = 
2G[a(n + A) + 2G] + A/t/3, a = s + 2Q and /3 = a 2 - 
4f2 2 . Here A = 57 — 1 is the atom-cavity detuning and 
(• • • ) denotes the asymptotic value. One may also obtain 
another useful expression for estimating the correlation 
between the atom and the field 

(Wi)=4-(w-w). (11) 

A + K 

For higher orders approximations the formulae are not 
so compact, so below we shall present only the numerical 
results for the 2-photons approximation. 

In many practical situations the circuit QED system 
contains other elements, apart from the selected cavity 
mode and the artificial atom. In particular, the atom al- 
ways interacts with other cavity modes, usually neglected 
due to their strong detuning from the atomic transition 
frequency and because initially they are not populated by 
photons. In ordinary cavities realized as coplanar waveg- 
uide resonators the modes are regularly spaced in fre- 
quency, and if the atom is fabricated at one end of the 
resonator, then all the cavity modes have a voltage antin- 
ode at the atom's position [?j]. To check whether other 
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cavity modes are relevant for the phenomenon of photon 
generation due to the ART, we suppose that, besides the 
selected cavity mode (the "true mode"), the atom also 
interacts with another cavity mode with frequency v and 
the corresponding coupling constant \fi>g - a "parasitic 
mode" , whose variables are denoted with tilde. In this 
case the total Hamiltonian becomes 0, [24| 

H = H c = Hr + vn + \fl>g pa y , (12) 

and assuming that the cavity quality factor is the same 
for the two modes, the corresponding modification in the 
Liouvillian is 

C c = C + vkD[o\. (13) 

Based on common experimental conditions, we assume 
that initially the parasitic mode is in the vacuum state, 
and consider two plausible scenarios, a) The true mode 
(v = 1) is the fundamental mode of the half-wave res- 
onator, and the parasitic mode corresponds to the first 
harmonic full wavelength resonance, v = 2. b) The oppo- 
site scenario: the true mode [y — 1) corresponds to the 
first harmonic full wavelength resonance, and the para- 
sitic mode is the half-wave fundamental mode, v = 1/2. 
Below we distinguish between these two cases by indicat- 
ing the value of v. 

Furthermore, in some circuit QED setups two almost 
identical artificial 2-level atoms are fabricated within the 
same cavity to perform two-atoms quantum gates, where 
the cavity mode plays the role of a quantum bus. When 
one is interested in performing only one-atom quantum 
gates, the second atom is effectively "turned off' by tun- 
ing its transition frequency fl far apart from the cavity 
mode's and the other atom's frequencies. So we study 
the scenario where the "connected" atom is the "true 
atom" and the "disconnected" one acts as the "parasitic 
atom" (its variables are denoted with tilde) , so the total 
Hamiltonian and the Liouvillian become [4fJ 

H = H A = H R + £lE + gpa y (14) 

C A = £+\-D[a-} + ^V[a z ], (15) 

where we assumed that the coupling constants and the 
dissipative rates are the same for the both atoms. We 
shall consider two scenarios: c) the parasitic atom tran- 
sition frequency is significantly smaller than the cavity 
mode frequency, = 0.2, or d) is significantly larger, 
h = 1.8. 

Below we study numerically the AMENDART of the 
true cavity mode, (n), and the true atom, (E), as func- 
tion of the system parameters under the 1- and 2- pho- 
ton approximations. We verified that in the weak cou- 
pling limit, 5 C 1, these quantities grow up approxi- 
mately quadratically in g, so in the figures [T] - 0] we set 
g = 5 x 10~ 2 , a value that can be achieved in some cir- 
cuit QED systems and is within the weak coupling limit. 



Moreover, as we are interested in the lower bounds for the 
mean excitation numbers, we set the reservoir tempera- 
tures to OK, n = 0. This simplification does not restrict 
the scope of our analysis, since we verified that for small 
but finite temperatures (n -C 1) the cavity's and atomic 
excitations due to finite temperature are approximately 
added to the AMENDART (data not shown). 

In the figure [1] we show how (n) and (E) depend on 
the true atom transition frequency fl in the presence of 
parasitic elements for parameters indicated in the cap- 
tion, that can be achieved in the near future experiments. 
The dashed lines correspond to the 1-photon approxi- 
mation and the solid lines - to the 2-photons approxi- 
mation. We see that for the chosen parameters the 1- 
photon approximation is quite accurate, and out of res- 
onance the AMENDART decrease as function of il. In 
the figure [2] we perform a similar analysis as function 
of the pure dephasing rate, from which we see that the 
AMENDART grow up as function of 7, and the 1-photon 
approximation loses its accuracy as (n) (or 7) increases. 
Remarkably, in these cases the parasitic elements can- 
not be neglected, since they contribute substantially to 
the AMENDART. Therefore, if one wants to know pre- 
cisely the values of AMENDART in circuit QED systems, 
one has to take into account the parasitic elements, even 
when they are far detuned and initiated in their respec- 
tive ground states. 
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FIG. 1: AMENDART as function of atomic transition fre- 
quency for parameters re = A = 1CF 6 and 7 = A/4 in the pres- 
ence of parasitic elements, whose parameters are indicated in 
the plots. Dashed lines correspond to the 1-photon approxi- 
mation, and the solid ones - to the 2-photons approximation. 

Despite a finite asymptotic photon population in the 
true cavity mode, it cannot be associated to an effective 
temperature reservoir (42j, since the resulting distribu- 
tion is significantly different from the thermal one. This 
is illustrated in the figure^, where we show the probabil- 
ity P n of having n photons in the cavity, obtained under 
the 2-photons approximation in the presence of the para- 
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FIG. 2: AMEND ART as function of 7 for parameters k = 
A = 1CT 6 , n = 1 in the presence of parasitic elements, whose 
parameters are indicated in the plots. The solid lines corre- 
spond to the 2-photons approximation, and the dashed ones 
- to the 1-photon approximation. 

sitic atom for parameters indicated in the caption (filled 
black bars). For comparison we also show the thermal 
distribution with the same mean photon number (red 
bars with sparse pattern), whereby one can see that the 
two distributions are quite different. Besides, there is a 
correlation between the "true field" and the "true atom" 
in the asymptotic state, as quantified by the quantum 
mutual information, that measures the total amount of 
correlation in a bipartite quantum state [431 ] 

I af = S(p a ) + S( Pf )-S( Paf ). (16) 

Here S (p k ) = — Tr[p fc log 2 (/Oj.)] is the von Neumann 
entropy of the kth subsystem whose dynamics is de- 
scribed by the reduced asymptotic density matrix p k with 
k = {atom, field, atom — field}. I a f is shown in the fig- 
ure [3Jd as function of f2 in the absence of the parasitic 
elements (black lines), and when the "parasitic atom" 
(orange lines) or "parasitic mode" (blue lines) are present 
for parameters K = A = 10~ 6 and 7 = A/4. The dashed 
(solid) lines correspond to the 1 (2)-photon approxima- 
tion, and the results for the 1-photon approximation in 
the absence of the parasitic elements can be obtained an- 
alytically with the aid of equations (f9)l - (p~T|) . One can see 
that asymptotically the atom-field system is correlated, 
Iaf i= 0, and the mutual information decreases in the 
presence of the parasitic elements, since in this case they 
acquire some information about the system of interest. 

One can also notice from the right plot of figure [3^ 
that when the true atom is out of resonance, (n) is larger 
for smaller k. This is better depicted in the figure |H 
where we show the total asymptotic mean excitation 
number under the 2-photons approximation, (712) + (E2), 



as function of damping rates under the resonance (sparse 
white curve) and out-of-resonance (filled colored curve) 
conditions in the presence of the parasitic atom, with 
the system parameters indicated in the caption. One 
can see that at the resonance (712) + (E2) does not de- 
pend strongly on the damping rates, contrary to out-of- 
resonance case, where it is smaller when n and A are 
approximately equal. 
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FIG. 3: a) Photon number distribution for parameters fl — 
0.2, A = 10~ 6 , 7 = A/4 and k and Q. as indicated in the 
plots. Filled black bars describe the photon distribution due 
to AMEND ART, and the red bars with sparse pattern cor- 
respond to the thermal distribution with the same (n). b) 
Quantum mutual information I a f for the asymptotic state 
with and without the parasitic elements. 



III. DISCUSSION AND GENERALIZATIONS 

The phenomenon of excitations generation above the 
thermal values occurs due to the combined action of the 
ART in the RH and the Markovian approximation for 
the reservoirs (36[. The ART describes continuous spon- 
taneous creation from vacuum and annihilation of vir- 
tual atomic and photonic excitations [23J, |24{ . Under the 
Markovian approximation, the reservoir correlation time 
is very short compared to the time scale for a significant 
change in the system [30| , i.e. the time needed to annihi- 
late excitations. Therefore, once the excitations are spon- 
taneously created, a fraction of them decays to the reser- 
voir due to atomic and/or cavity damping before they can 
be annihilated, so the system coherence is destroyed and 
the initial zero-excitation state cannot be restored. From 
the point of view of quantum trajectories [H, [3(J IHj] , the 
continuous monitoring of the system by the Markovian 
environment promotes the virtual excitation into real 
ones due to the continuous change of system's state, while 
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FIG. 4: Logarithm of the total mean excitation number 
log((n2) + {E2)) as function of log(ft) and log(A) in the pres- 
ence of the parasitic atom, for parameters 7 = A/4 and 
0, = 0.2. The colored curve corresponds to out-of-resonance 
case, = 0.7, and the white one to the resonance, fi = 1. 



a fraction of the excitations is destroyed due to destruc- 
tive measurements, thereby asymptotically nonzero pho- 
ton and atomic populations are built in the system. The 
pure atomic dephasing increases the number of created 
excitations (37j because the dephasing reservoir performs 
nondestructive measurement of the atomic state [30| , so 
the virtual excitations are promoted into real ones nonde- 
structively. Besides, the Markovian pure atomic dephas- 
ing can be attributed to random changes of the atomic 
transition frequency [H, [33l , which promote the vir- 
tual excitations into real ones due to the modification of 
the system ground state, analogously to the dynamical 
Casimir effect p^l - fr7l . |44| . It is worth noting that, by def- 
inition, the pure dephasing reservoir always has a finite 
temperature [3fj| , whereby it can perform nondestructive 
measurements of the atomic state, so the substantial am- 
plification of system excitations due to pure Markovian 
dephasing is not so puzzling. 

The phenomenon of nonzero asymptotic photon pop- 
ulation also occurs for a single dissipative cavity if one 
takes into account the anti-rotating cavity-reservoir in- 
teraction. Suppose we isolate a single constituent of the 
reservoir (e.g. 2-level atom, harmonic oscillator, etc) and 
treat it as an ancilla; then, one performs the standard 
Born-Markovian master equation treatment over the re- 
maining reservoir constituents. Thereby we obtain a 
Markovian master equation for the cavity-ancilla system, 
including the anti-rotating interaction. By performing 
the asymptotic analysis of the previous section and trac- 
ing over the ancilla variables at the end, one would end 
up with a nonzero photon population due to the virtual 
cavity-ancilla excitations being promoted into real ones 
by the Markovian reservoir. Another way of arriving at 



this conclusion is to use the simplified phenomenological 
approach. The mo st g eneral master equation for a single 
cavity field mode |45l l46j . preserving the normalization 
and hcrmiticity of the statistical operator p and contain- 
ing only bilinear forms of operators x and p is given by 
the equation © with the effective Hamiltonian 



H = H 



(17) 



where Hq = n is the cavity Hamiltonian (recalling that 
the cavity frequency is set to v = 1), and the damping 
superoperator is C = C' K 



£'kP = -j ({P,{x,p}}- [x,{p,p}}) - D p [x,[x,p}} 

- D x \p, \p, p}} + D z {[x, [p, p]} + [p, [x, p]\) , (18) 

with p, k, D x , D p and D z being arbitrary time- 
independent coefficients under the Markovian approxi- 
mation. The condition 



D p D x 



Dl > ( K /4) : 



(19) 



guarantees that the positivity of the statistical operator 
is preserved for all times and for any physically admis- 
sible initial state. This is the necessary and sufficient 
condition (together with conditions D x > and D p > 0) 
of rcducibility of the superoperator (fT8"|) to the Lindblad 
form (4?| ■ The standard master equation (|3]) corresponds 
to the choice D p — D x = k (1 + 2n) /4 and D z = p = 0, 
which satisfies the inequality (|19l) . 

One can verify (36| that asymptotically the vacuum 
state (with mean values (x 2 ) = (p 2 ) = 1/2, ({x,p}) = 0) 
is achieved for the coefficients 



D' p = {K + p)/A, D' x = ( K -p)/A, 



D' z - 0. (20) 



Using the inequality (IT^|) one obtains n 2 — p 2 > k 2 , whose 
solution is p, = 0. Therefore the vacuum state can be 
achieved only for the coefficients D' x = D' p = k/4 and 
p = D' z = 0, which correspond to the standard master 
equation © at zero temperature. However, the SME 
is deduced microscopically by making the RWA on the 
cavity-reservoir interaction [23l. l24l 130] . i.e. by neglect- 
ing the anti-rotating terms responsible for simultaneous 
creation of one virtual photon and a virtual reservoir ex- 
citation. Therefore, if the anti-rotating cavity-reservoir 
interactions are taken into account, under the Marko- 
vian approximation the asymptotic mean photon num- 
ber in the dissipative cavity is larger than zero. Finally, 
if one uses a single dissipative kernel (fT8"|) together with 
the Hamiltonian (|T7j) to describe the circuit QED sys- 
tem, with Hq = Hr, under the 1-photon approximation 
one gets (for g ^ 0) 



(Ei) 



2 I 2D 



■i 1 1 1 2(j 2 D x 

1 nQ,(j)v + D X p/ Dp 

2 D xp + 2giD x( f> 



> 



>o, 



(21) 
(22) 
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where D xp — D x + D p , v± = (1 ± 2D Z ) and (f> = 

[O 2 + 4£> 2 + v + v-D x / 'Dp] . This demonstrates that 
regardless of the exact form of the master equation and 
the amount of dissipative channels, under the Marko- 
vian approximation the asymptotic mean photon and the 
atomic excitation numbers are always greater than zero, 
so the vacuum state is never achieved exactly. 

IV. SUMMARY 

In summary, we studied the behavior of the cavity 
field's and atomic asymptotic mean excitation numbers 
due to anti-rotating term (AMENDART) for typical pa- 
rameters of the Markovian circuit QED system, showing 
that "parasitic elements" (such as other cavity modes and 
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